ON THE EFFECTIVE LIGHT-CONE QCD-HAMILTONIAN 



H.C. PAULI 

Max- Planck- Institut fur Kernphystk, D-69029 Heidelberg, Germany. 
E-mail:pauli@mpi-hd. mpg.de 

Taking the effective interaction between a quark and an anti-quark from previ- 
ous work, the dependendence on a regularization scale is removed in line with the 
renormalization group. In order to emphasize the essential point, the full spinor 
interaction is replaced by a model which includes only the Coulomb and the hyper- 
fine interaction. By adjusting the effective quark masses, the only free parameters 
of the theory, the mass and the size of the pion are reproduced, as well as the mass 
of all other pseudo-scalar mesons. Estimates for the vector mesons are close to 
the empirical values. The model exposes screening rather than strict confinement. 
The ionization thresholds are in general much larger than the pion mass. 



1 Introduction 

The bound-state problem in a field theory particularly gauge field theory is 
a very difficult many-body problem and has not been solved thus far. But 
perhaps an exact solution is not necessary for reconciling the so succesful con- 
stituent quark models with a covariant theory such as QCD and for being 
able to compute observables like the structure functions and distribution am- 
plitudes from an underlying theory. Perhaps it suffices to have approximate 
or QCD-inspired solutions which are not wrong from the outset. The present 
work is still ongoing and contributes to these questions. Perhaps one begins 
to understand the essential point. 

2 A QCD-inspired effective interaction 

The light-cone approach to the bound state problem aims at diagonalizing the 
invariant mass-squared operatorlll Its matrix elements are obtained directly 
from the .gauge field Lagrangian in the light-cone gauge and are explicitly 
tabulated!!! The many-body aspect of the problem was reduced recently, by 
the method of iterated resolvents p to the problem of phrasing an effective 
Hamiltonian which acts only in the Fock space of one quark and one anti-quark. 
The reduction is not exact but rests on certain simplifying assumptions which 
can be checked only a posteriori. Quite cautiously, one should therefore speak 
of a 'QCD-inspired' effective Hamiltonian. Important is that both the effective 
and the full Hamiltonian have in principle the same eigenvalue spectrum and 
that the Fock space is systematically reduced to the gg-space. Important is 
also that the Fock-space is not truncated, that all Lagrangian symmetries 
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Figure 1: The three graphs of the effective 
interaction in the qq-space according top. The 
upper graph illustrates an effective interaction 
due to a two gluon annihilation. For flavor-off- 
diagonal mesons its amplitude is kinematically 
suppressed. The lower two graphs correspond 
to an effective one gluon exchange. — The effec- 
tive interaction scatters a quark with four mo- 
mentum and helicity (fci, Ai) into (fcj, A^), and 
correspondingly the anti-quark from (^21-^2) 
into (k' 2 , Aj). 



are preserved, and that the higher Fock-space amplitudes can be retrieved 
systematically from the qq-projection once that is available. 

The effective interaction has several contributions which are illustrated 
in Fig. |l|. The present work is restricted to ffavor-off-diagonal mesons, to 
mesons where quark and anti-quark have different flavors. Then the effective 
interaction due to the annihilation of two gluons is kinematically suppressed, 
and the final (one-body) integral equation in the qq-space is governed by an 
effective one gluon exchange, i.e. 
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Here, M 2 is the invariant-mass squared eigenvalue and (x, k±; Ai, A2] V 1 ) t ne 
associated eigenfunction. It is the probability amplitude for finding a qq-Fock 
state in which the quark has longitudinal momentum fraction x and transversal 
momentum k± and the anti-quark correspondingly f — x and — and where 
the respective helicities are Ai and A 2 . The mean four- momentum transfers of 
the quarks and the spinor factor are respectively defined by 



= -\[(ki-k[) 2 + (k 2 -k / 2 ) 2 ], 



(2) 



S = (A 1 ,A 2 |5|A' 1 ,A 2 ) = [u(fc 1 ,A 1 ) 7 ^(^,A' 1 )] [v(k' 2 , X 2 ) lfl v(k 2 , X 2 )} . (3) 

The effective quark masses mi j2 and the effective coupling constant a(Q; A) 
depend both, in general, on a regularization scale A. The regulator function 
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i?(x',/c^; A) restricts the range of integration and in general is a function of 
the same A, see also below. 

The current work is restricted to the lowest order of approximation (LOA) 
where m — m and a(Q; A) = ql The expressions for the next-to-lowest order 
(NLO) can be found elsewhereo One has a fair confidence jnto the validity of 
Eq. ([!]), since the alternative method of Hamiltonian flowB leads to the very 
same eauation as in LOAq The same equation had also been obtained priox 
to thatpEl with however less stringent arguments. The work of Wilson et al. 
had the same emphasis, but that it did not lead to the same formulas, see also 



3 The crucial point in a simple model 

Why is a regularization necessary at all? 

In light-cone paramctrization, the quark is at rest relative to the anti- 
quark when k± = and x — x = mi /(mi + m^). For very small deviations 
from these 'equilibrium values' the spinor factor is diagonal in the helicities, 
(Ai, A2 1 I A 1 , X' 2 ) ~ 4toiTO2 Sx 1> x / 1 S\ 2 ,y 2 i as can t>e veryfied from an explicit 
presentation of S% For very large deviations, particularly for k^ 2 ^> k 2 , holds 
Q 2 ~ k 12 and (tl |5| |4) — 2fc{ 2 . The ratio is then a simple dimensionless 
number, S/Q 2 — 2, independent of fc^. Both extremes are combined here by 

(U|5|U) 2 _ 4mim 2 

Q ~ Q 2 [) 

which is considered as a model of the interaction in singlet states, dropping the 
less important spin-orbit interaction. nn 

Unfortunately, I have not realized earlier Era that the innocent and hnitc 
number '2' in Eq.(||) has the same orign as the familiar divergencies in gauge 
field theory. The latter arise, typically, when bilinar products of the elemen- 
tatry Dirac interaction like in S divided by some energy denominator like Q 2 
are integrated over all phase space. 

It was ijlear from the outset that the many-body Hamiltonian has to be 
regularizedEl Since fancy methods like dimensional regularizatioix are not ap- 
plicable, Fock-space regularization was considered to be sufficientEl In practice, 
it has lead to difficulties and was was replaced by vertex regularizationp where 
the elementary Dirac interaction at a vertex was endorsed with kind of a form 
factor F(A), i.e. 

[u(l)^u(2)} E/ ,(3) — [«(l)y«(2)] e„(3) F(k u k 2 , k 3 ; A). (5) 
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Requiring that the off-shell mass of the two scattered particles Mq = (k 2 + k^) 2 
is limited, F becomes the step function F(k\, hi, ks', A) = Q(M§ — A 2 ). In 
the absolute squares of the effective interaction, as in Eq.(Q), it appears as a 
regulator function R as a consequence of regulating the theory from the outset. 

4 Rewriting the integral equation in conventional momenta 

Unpleasant is that x and k± have different ranges (0 < x < 1, — oo < k x ^.oo). 
Instead of x a k z ( — oo < k z < oo) can be introduced by the transform!!!] 

Ex + k z 



X ( k *) = TP -L TP 1 With E ^ = -2 + k l + ^ ( 6 ) 

The Jacobian of the transformation can be cast into the form 

dx 1 dk z , . 1 EiE 2 , , 

with A(k) = — - \ I , (7) 



x(l — x) A(k) m r ' to,. Ei + E2 

and = — — h = -, m s = mi + TO2 ■ (8) 

m r TOl TO2 

In the same convention the kinetic energy becomes 

= mj + kl + to 2 + kl _ ^ ^ p 
a; 1 — a; 

with C(fe) = (£1 + mi + £ 2 + ffh) [ -= * + ) , (10) 

+ toi E 2 + m 2 J 

without approximation. If one substitutes the wave function according to 



^A) = J^|0(iA) (ii) 

y a;(l — xj 

and drops explicit reference to the helicities, one arrives at an equation, 

M2-mS-C(fc)P]^V- * / 7 £P= (12) 

TO r J . I A(U\ AIU\ <* 



'A(jfe)A(jfc 

which is identical with Eq. (Q) , but which looks like one with usual 3- momenta. 
5 Interpretation in configuration space 

The main reason for introducing conventional 3-momenta is interpretation. It 
is more transparent in configuration space. Since the various factors A{k) 
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with their square-roots prevent straight-forward Fourier transforms, the non- 
relativistic approximation (fc' 2 <C m\ 2 ) is applied consistently. Using 



A(k) = l, C(k) = —, Q 2 = (k-k') 2 



(13) 



the model interaction of Eq.(|J), and in addition the regulator R = 1, Eq.(^2|) 
leads to a local Schrodinger equation 



MZ 2 ""Sri 
— m s k 

M 2 -m 2 s - 
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<KV), (14) 



?/>(r) = 2m s V(r)ip( , r), with p = iV. 



(15) 



The potential V(r) bears great similarity with 14 1 £_ g (r), the hyperfine interac- 
tion in the singlet channel found for hydrogen in all textbooks, i.e. 



V(r) = ~a(- + 2 



— ^>(r)V 
i r m s J 



V M _ s (r) = - a(-+ g p S^(r) 

4 r m e m p 



(16) 
(17) 



The strange '2' in Eq.(|4|) finds its explanation as the gyromagnetic ratio for a 
fermion, with g p = 2. But Eq.(p^) has no solution! A Dirac-delta function is 
no proper function and must be regulated, for instance by a Yukawa-potential 
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(18) 



Important is that the delta (J d 3 fS^(r) = 1) and the Yukawa have the same 
strength J d 3 f |/t 2 exp(— /xr)/(27rr)] = 1. Transforming back to momentum 
space gives 



l\ r 2 2 "'s ,2 

M — m. k 



^ id 3 k> f ^ 



3tt 2 



1 



Q 2 2m r n 2 + Q 5 



0(fc')-(19) 



Obviously, replacing in configuration space the Dirac-delta function with a 
Yukawa potential corresponds to introducing in momentum space the regulator 
function R = [i 2 j{\i 2 + Q 2 ). 

The regularization of a Dirac-delta .function is an old theme of nuclear 
physics in the context of pairing theorytil It was the.-Boint of orign for the 
similarity transform( 12 H 13 l and was investigated recentlylp again. 



Figure 2: Coupling constant a versus reg- 
ularization scale fj, in units of 350 MeV. 



Figure 3: The spectrum of the first three 
singlet-s states in MeV versus regulariza- 
tion scale /i in units of 350 MeV. 



6 Renormalization 

The eigen values of the integral equation Eq. ( |l9|) depend now on a regulariza- 
tion scale fi. In line with the current understanding of field theory one replaces 
a and m by functions of /u, such that every eigenvalue is stationary against 
small variations of a, m and fi, i.e. dM^/dfi — 0! The general procedure, 
however, is not possible here, since a and fi appear in Eq.(|l9|) in the typical 
combination a/j, 2 /(fi 2 + Q 2 ). One therefore proceeds by looking for a function 
such, that the calculated mass of the pion, for example, agrees with the 
empirical value. The so obtained function is considered universal. The 
actual value of fj, can then be fixed by a second requirement. 

7 The pion 

For carrying out, in practice, the programme of Sec. [(] one needs_an efficient 
tool for solving Eq.(|l9|). Such one has been developed recentlyfj restricted 
however for simplicity to spherically symmetric s-states. Fixing the up and 
the down mass to rn u = Wid = 1.16 (the unit of mass is chosen here as 350 
MeV/c?), the calculated pion mass is required to agree with the experimental 
value 113 to within eight digits. For every value of fi one obtains thus an 
which is displayed in Fig. |^. The even greater surprise came when investigating 
the spectrum of the pion along the so obtained a M , as given in Fig. |. The 
lowest state, the 7T + , stays nailed fixed to the empirical value, while the second 
and the third (as well as the higher ones) have an extremum at [i ~ 3.8. Such 
an extremum was not really expected, but in view of the renormalization group 
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Figure 4: The pion wave function < I > (p) 
is plotted versus p/ 1.552 (filled circles) in 
an arbitrary normalization. The diamonds 
denote a Coulomb wave function with the 
same a and the open circles an approxi- 
mate wave function, see text. 
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Figure 5: The relativistic potential energy 
W(r) for the pion is plotted in GeV ver- 
sus the radius in fm (solid line). The pion 
mass is indicated by the dashed, and the 
ionization threshold by the dashed-dotted 
line. 



is useful: /j, determines itsself from the solution! 

For sufRently large values of [i, say for fi 3> 10, see Fig. ||, the coupling 
constant decreases strongly and the spectrum becomes more and more the 
familiar Bohr spectrum with a small hyperfine shift. Obviously, it can be 
calculated perturbatively, indeed, for sufficiently small a. 

Having fixed now the regularization scale fi ~ /xrj = 3.8 gives the coupling 
constant a ~ a,, = 0.6904, see also Fig. |^. The wave function of the pion is 
plotted in Fig.Hl By a fit to an approximate expression <& a (p) = 10/(H-p 2 /p 2 ), 
a width p a = 1.471 is obtained which is much larger than p c = |am r = 0.5339, 
the Bohr momentum of the Coulomb wave function with the same a. Having 
the wave function, one could, in principle, calculate theJorm factor and the 
distribution amplitude according to the exact expressionscl This is however not 
a trivial undertaking, and, instead of, the above approximate expression was 
Fourier transformed to ip(r) ~ e~ PaT . The calculated root-mean-sqare radius 
of the pion becomes then 

< r 2 >^= — = 0.663 fm. (20) 

Pa 

The good agreement with the experimental value < r 2 >i = 0.657 fh£3 is 
actually the result of fitting the quark masses m u = m u — 406 MeV, see 
above. This exhausts all freedom in choosing the parameters of Eq.([l9|). 
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One should emphasize that the Yukawa potential in Eq. (|19|) at the scale /io 
pulls down essentially only one state. The others are left rather unperturbed 
at their Bohr values, see Fig. [| The_analogy with the pairing model in early 
nuclear physics is more than obviouslij In a future and more rigorous solution 
of the full Eq.([l2]) one expects that the excited states can be disentangled into 
almost degenerate singlets and triplets, which in turn can be interpreted as an 
excitation of the pseudo-scalar pion, or the ground state of a vector meson, 
respectively. In any case, the first excited state of the simple model correlates 
very well with the mass of the p + and the other vector mesons, see below. 

8 On the question of confinement 

The question of confinement is closely related to the potential V(r), but the 
relation is subtle. One of the great advantages of light-cone quantization is the 
additivity of the free part and the interaction. Because of the dimensions of 
invariant mass squares the relation of potential and kinetic energy is somewhat 
hidden. In analogy to a classical Hamiltonian, the invariant mass-squared op- 
erator M 2 in Eq. ( |19| ) is interpreted as the square of a position and momentum 
dependent mass 

M(r,p)=m s Jl + ^— + —V(r). (21) 
W m s m r m s 

For V(r) <C m s and p 2 <C m s m r it could be expanded as 

2 

M(r,p) ~ m s + ^— + V(r), 

~ rest mass + kinetic energy + potential energy, 
These conditions, however, can not be satisfied everywhere. I therefore propose 
to introduce a relativistic potential energy by W(r) = M(r,p — 0). This 
function is plotted in Fig. ||. It vanishes at vq = 0.5194 fm and the classical 
turning point is r t — 0.5300 fm, in agreement with Eq. (pcf) . Since the sum 
of the quark masses is 812 MeV, the ionization threshold occurs at 542 MeV: 
In order to liberate the quarks in the pion, one has to invest more than three 
times the rest energy of a pion. 

9 The meson masses 

The remaining parameters of the theory, the masses of the strange, charm 
and bottom quark can then be determined by reproducing the masses of the 
pseudo-scalar mesons K~ D° and B~ respectively. This gives fn s = 508, 
m c = 1666 and nib = 5054 MeV. The remaining off-diagonal pseudo-scalar 
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Table 1: Left: Experimental masses of the flavor-off-diagonal physical mesons in MeV. 
Scalar mesons are given in the lower, vector mesons in the upper triangle. - Right: The 
calculated mass eigenstates of QCD in MeV. Singlet- Is states are given in the lower, singlet- 
2s states in the upper triangle. 



mesons are then calculated straightforwardly and compiled in Table |lj. In 
view of the simplicity of the model, the agreement with the empirical values is 
remarkable, indeed, and at least as good as for any other model. The strong 
correlation of the first excited state with the vector mesons was mentioned. 

10 Summary and some conclusions 

The very difficult problem of solving the many-body Hamiltonian for QCD is 
replaced here by the problem of solving an effective QCD-inspired Hamiltonian. 
It is obtained from the QCD-Lagrangian in the light-cone gauge by the method 
of iterated resolvents or the Hamiltonian flow equations in the lowest non- 
trivial approximation. In order to avoid the numerical complexities of the full 
solution in light-cone coordinates, this effective interaction was stripped-off of 
almost all ingredients except the Coulomb and the hyperfine interaction, which 
are analyzed in a non-relativistic treatment with only 3-momcnta. It is argued 
why the hyperfine interaction in the singlet channel needs to be regulated like 
the usual divergencies in gauge field theory. 

It was courageous to apply such a simple model to the mystery particle 
of QCD, to the pion, but it turns out that both the pion's mass and size 
are obtained without vacuum condensates of any kind, just by adjusting the 
canonical parameters of a gauge field theory, the coupling constant and the 
quark masses. The necessary regularization constant determines itself by the 
theory, in line with the current interpretation of the renormalization group. For 
the pion, as well as for the other mesons, one gets thus for the first time a QCD- 
inspired light-cone wave function by which one can compute the form factor 
by the available exact formulas, as well as the higher Fock space amplitudes. 

The present approach exposes screening rather than confinement. The 
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removal energy of a quark from a meson is higher than the pion's mass. A bare 
quark is thus unlikely to be observed. One wonders whether this is not much 
more a physical picture than the strict linear confinement obtained by lattice 
gauge theory, the only available alternative approach for a non-perturbative 
description. But like the light-cone approach, the latter is burdened with its 
own pecularities among them the far extrapolation down to the pion. 

In fact the present model runs short in several aspects. That the regulator 
function should apply to the whole kernel, and not only to one piece of the 
interaction, is a technical detail, which can be easily corrected in future work. 
More importent is whether the two-gluon annihilation interaction, which was 
omitted here on perpose, can account for the flavor diagonal mesons and the 
so important aspects of isospin. 

Astounding is also that the present model for QCD differs from QED 
only by the color factor | in the coupling constant. Can it be true that the 
essential difference between these two theories with their so distinctly different 
phenomenology is manifest only as large and small coupling? 

Depite all of that I conclude that the light-cone approach to hadronic 
physics may now have a fair chance to compute the hadronic spectra and to 
predict the hadronic wave functions in line with a covariant theory. 
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